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The integrability of the one- dimensional (ID) fermion chain model is investigated in the 
framework of the Quantum Inverse Scattering Method (QISM). We introduce a new R- 
operator for the fermion chain model, which is expressed in terms of the fermion operators. 
The i?-operator satisfies a new type of the Yang-Baxter relation with fermionic L-operator. 
We derive the fermionic Sutherland equation from the relation, which is equivalent to the 
fermionic Lax equation. It also provides a mathematical foundation of the boost operator 
approach for the fermion model. In fact, we obtain some higher conserved quantities of 
the fermion model using the boost operator. 
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1 Introduction 

In the last decades, many integrable spin chain models have been investigated by means of 
Quantum Inverse Scattering Method (QISM for brevity). |l| In the QISM, the i?-matrix, 
which satisfies the Yang-Baxter relation with the L-operator 

12 2 1 

R 12 (u 1 ,U2)Lj(ui)Lj(u 2 ) = L j (u 2 )L j (ui)R 1 2(u 1 ,u 2 ), (1) 

plays the essential role, [fj, [| f|, ^ ^[ For the spin chain models, the i?-matrix Ri 2 (ui, u 2 ) 
is a c-number matrix and satisfies the Yang-Baxter equation 

R 12 (u 1 ,u 2 )R 13 (u 1 ,u 3 )R 23 (u 2 , u 3 ) = R 23 {u 2 , u 3 )R 13 (ui, u^R^ux, u 2 ). (2) 
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From the Yang-Baxter relation ([I]), it follows that the transfer matrix defined by 



t(u) — ti a (L N (u) . ..Li(u)) (3) 
constitutes a commuting family 

[t(u),t(v)}=0. (4) 

Expansion of t(u) in terms of spectral parameter u gives the conserved quantities. There- 
fore, this proves the existence of the conserved operators which are involutive, 

[/CfljW] = . (5) 

We now know that there also exist many solvable fermion chain models which are 
important in condensed matter physics. They are usually related to the spin models 
through the Jordan- Wigner transformation. For example, the fermion chain model related 
to the XX Z model is given by |7|, |8j 

N ^ 

% = zZ{ c h^ + 4+i c ^ + y( 2n i - 1 )( 2 ^+i " !)}• ( 6 ) 

3=1 

Here cj and Cj are the creation and the annihilation operators at site j which satisfy the 
canonical anti-commutation relations 

{cj, c k } = {ct, c[} = 0, {cj, c[} = 8 jk , (7) 

and rij is the number operator 

Tlj = c]cj. (8) 
The periodic boundary condition (PBC) is assumed in (|J), 

Cn+1 — Ci, C N+1 = C l- (9) 

We refer to this model (§) as the XX Z fermion model. 

Since the Jordan- Wigner transformation connects the spin model and the fermion 
model, physical properties of the fermionic chain models including the exact integrabil- 
ity are often discussed using the corresponding spin chain model. However, there are some 
important differences between the fermion models and the spin models. In particular, the 
PBC for the fermion model does not correspond to the PBC for the spin model due to 
the non-locality of the Jordan- Wigner transformation. Therefore it is more appropriate to 
treat the fermion models keeping their fermionic nature as much as possible (c.f. refs. 9 
and 10). 

The QISM for the fermion chain model (||) has been developed by several authors, 
|8], [IT], where the fermionic L-operator is introduced. The fermionic i?-matrix which 



2 



satisfies the graded Yang-Baxter relation was also found. The existence of the fermionic 
i?-matrix can be used to establish the integrability of the fermion model under the PBC (§) 
or other twisted boundary conditions. The most intriguing point of this method is that the 
matrix elements of the fermionic L-operator consist of the fermion operators, which means 
that the quantum space is fermionic. On the other hand, the auxiliary space remains to 
be a usual spin space and correspondingly the f?-matrix has the c-number elements. The 
fact suggests that there still requires further investigation on the fermionic formulation of 
the QISM. 

In this paper, we introduce an i?-operator for the XXZ fermion model. It consists 
of the fermion operators and satisfies a new type of the Yang-Baxter relation with the 
fermionic L-operator. The role of the auxiliary space and the quantum space for the 
fermionic L-operator is exchanged in the Yang-Baxter relation. We derive the fermionic 
Sutherland equation from the Yang-Baxter relation, which is shown to be equivalent to the 
fermionic Lax equation. O, [TB| To obtain the higher conserved operators for the lattice 



spin models, it is in general very useful to introduce the boost operator. |L4], [L5|, [16] The 



fermionic Sutherland equation provides the basis for the application of the boost operator 
to the fermion models. We derive some higher conserved operators for the XXZ fermion 
model using the boost operator. 

The paper is organized as follows. In §2, we summarize salient points of the QISM for 
the XXZ fermion model. In §3, we consider the exchange of the auxiliary and quantum 
spaces of the fermionic L-operator and introduce an i?-operator. We also discuss the 
fundamental properties of the i?-operator. In particular, we find that the /^-operator 
satisfies the Yang-Baxter relation. In §4, we further discuss some applications of the 
fermionic L!-operator and the Yang-Baxter relation. We derive the Sutherland equation 
and introduce the boost operator. The last section is devoted to concluding remarks. 



2 Graded Yang- Baxter Relation for the XXZ Fermion 
Model 

In this section we briefly summarize the QISM for the XXZ fermion model (^|). The 
fermionic L-operator is given by 

jK ) V -2i/3(«)ct a(u)(l- nj ) -yyi^nj ' K ' 



where 



a{u) = sin(w + 2rj), 
7(u) = sin it, 

2(3{u) = sin2??. (11) 
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We express by (g> the Grassmann (graded) direct product 

s 

[A®B] a ^ 5 = {-ir^ p ^)A aP B, s 

s 

P(1) = 0, P(2) = l. 



(12) 



Then there exists the fermionic P-matrix which satisfies the graded Yang-Baxter relation 
(Fig. 1), 



12 2 1 

fti 2 (wi - M 2 )£(ui)£(m 2 ) = L(-u 2 )L(-ui)^i 2 (-ui - u 2 ) 
where, with J being 2x2 identity matrix, 

1 2 

= iCj(wi) ® J, £(w 2 ) = J <8> £jj(u 2 ). 

s s 

The explicit form of the fermionic P-matrix is 



(13) 



(14) 



ftia(tt) 



/sin(w + 2^) 

— isin-u sin 2?? 

sin2^ isin-u 

V o oo 



(15) 



— sin(u + 2r/) / 
It fulfills the graded Yang-Baxter equation (Fig. 2) 

ft 12 0)ft 13 (« + w)ft 23 0) = ^ 23 (^)^i 3 (n + u)tt 12 (u) 
Note that non-zero elements of ( |i~5|) are even with respect to the parity P(a), i.e., 
P(a) + P{(3) + P(a') + P((3') = (mod 2) for ft a/W /3' ^ 0. 



(16) 



(17) 



The monodromy matrix T(u) is defined as an ordered product of the fermionic L- 
operators 

T(u) = £#(«)... £i(u). (18) 
From the local relation (|13"D , we have the global relation for the monodromy matrix 

(19) 



12 2 1 

Ru(ui - u 2 )T j (u 1 )T j (u 2 ) = Tj{u 2 )T j {ui)'R 1 2{ui - u 2 ). 



The transfer matrix r(u) is defined by the supertrace of the monodromy matrix T(u) 

t(u) = strT(V) = ti{a z T(u)}. (20) 

Then the global graded Yang-Baxter relation (|19"D leads to the commutativity of the transfer 
matrix, 



[t( Ui ),t(u 2 )]=0. 



(21) 



The last equation establishes the exact integrability of the model (j^) under the PBC 
fact, we have 

t(u) = t(0){1 + u% + ...}, 

with the identification 

A = cos 2rj. 

For later use, we introduce a notation 

N 

where 

Kj,k = ^~^{( c i Cfc + c k c i) + \ cos 277(2^- - l)(2n fc - 1) + ^ cos 277}. 
We call IKj^ the two-point Hamiltonian density. 





Fig. 2 Graded Yang-Baxter equation, eq. (16) 



3 Fermionic R- Operator 



One of the characteristic features of the fermionic L-operator fllTf ) is that its matrix ele- 
ments consist of the fermion operators, that is, quantum space of the fermionic L-operator 
is fermion Fock space. However, the auxiliary space remains to be a usual spin space. The 
fermionic L-operator is therefore asymmetric with respect to the exchange of the quantum 
space and the auxiliary space. The reader may recall the similar situation when we con- 
struct the spin-S 1 XXZ model from the spin-| XXZ model. In this case, we first consider 
the intermediate L-operator, whose quantum space consists of spin operators of magnitude 
5*. Then we exchange the auxiliary and quantum spaces of the L-operator and consider a 
new Yang-Baxter relation. The i?-matrix which satisfies the Yang-Baxter relation should 
be in a (25 + 1) x (25* + 1) matrix form. From the f?-matrix, we can construct a new 



transfer matrix which generates the spin-S* XXZ model. || [17], 

We apply the similar procedure to the fermion chain model. Consider the following 
new Yang-Baxter relation (Fig. 3) 

a a a a 

ft f 12 («i - u 2 )Li(u 1 )L 2 (u2) = £ 2 (w2)£i(wi)ft f i 2 ( M i - u 2 ) (26) 

Notice that the role of the auxiliary and quantum spaces of the fermionic L-operator is 
exchanged in (|T^) and fl2"6]). Inevitably, the i?-operator Jl f (u) is no longer c-number matrix 
and consists of fermion operators. Of course, there is no assurance for the existence of such 
i?-operator at this stage. We assume a form of the i?-operator as 

0i[ 2 (ui - u 2 ) =gin 1 n 2 + g 2 n 1 (l - n 2 ) 

+ #3(1 - ni)n 2 + #4(1 - ni)(l - n 2 ) 

+ 95c\c 2 + g&c\c u (27) 



where g i} i = 1, . . . ,6 are scalar functions of ui and u 2 . Substituting the expression ( P7|) 
in (|26|), we have the following equations, 

94 = -9i, 93 = 92, 96 = 95, 

oc{u 2 )g x = -iiu 2 )g 2 - a(Mi)^ 5 , 
l(ui)9i = -a(u 1 )g 2 - 7(^2)^5, 
(3(u 1 )(3(u 2 )g 2 = {j(ui)a(u 2 ) - a(ui)^(u 2 )}g 5 . (28) 



A set of functional equations (|28| ) is solved to give (up to an overall factor) 

sin(ni — u 2 + 2r)) 



91 = ~9a 

92 = 93 = 



sin 2rj 
sin(-u 1 — u 2 ) 



sm2r] 

95 = ge = l. (29) 
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Thus obtained i?-operator, which we call fermionic /^-operator, 



\K n {ui - u 2 ) = — -{1 - Tlx — n 2 ) 



sin 2r] 
sin(wi — u 2 ) 
sin 2 77 

+ c\c 2 + 4ci, 



+ 



(ni + n 2 - 2n x n 2 ) 



satisfies 

The fermionic i?-operator enjoys the following properties. 
1. Regularity (Initial condition) 



%\ 2 {U = 0) = %X2- 



Here 



0C 



l-( C |_ C t)( C ._ C .) 



is the permutation operator for the fermion operators: 



DC 



J.I 



0C{jCi CjOCijj OCijC^ CjOCij, 



or or 



1. 



2. Local Hamiltonian: 



dX[ 2 (u) 



du 



u=0 



0Cx 2 0~C 



12- 



Here 'Ki 2 is the 2-point Hamiltonian density given in (pop. 
3. Yang-Baxter equation: 

The fermionic i?-operator satisfies the Yang-Baxter equation (Fig. 4) 



# f l 2 K - «2)ft f l3(Wl)ft 23 M = ^M^Wl^K - U 2) 



4. Unitarity: 



sin 2 u 



sin 2r). 



1. 



Now we consider a product of the fermionic i?-operators 

7l(u) = Xi N (u)...Xl 2 (u)3l { al (u). 
Then from (|35|), we see that 7^ (it) satisfies the Yang-Baxter relation, 

^ f 12 («! - « a )^(ui)a|(ti 2 ) = T\{u 2 )T{{ui)&[ 2 {ux - u 2 ). 



Here we note the commutativity, 

3^K)34M = 34M^>i), U k). 



(39) 



We refer to 7j{u) as monodromy operator. We can also introduce an analog of the transfer 
matrix, 



t\u) = StrO^(u) 

= a <o|ai(u)|o> - a (i|ai(«)|i> 8 , 



(40) 



where |0) a and |l) a are defined by 

C a \0) a = 0, |l) a = C t|0) a . (41) 

The fermionic transfer operator t { (u) constitutes a commuting family due to (p8[) 

[r i (u),r i (v)} = 0, (42) 

which also gives Hamiltonian ([]), 

N n J / 



ar_X^nr _ d_bgT> 



w=0 



(43) 





Fig. 3 Fermionic Yang-Baxter relation, eq. (26) 





Fig. 4 Fermionic Yang-Baxter equation, eq. (35) 
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4 Sutherland Equation and Boost Operator 

In this section we discuss further some applications of the fermionic i?-operator and the 
Yang-Baxter relation. 

First, we differentiate the Yang-Baxter relation (p6| ) with respect to the spectral pa- 
rameter u 2 at u 2 = u\. Then we get 

-% l2 'K 12 L l {u 1 )L 2 {u l ) + % 12 L l {u 1 ^ LA "- ) 



du 2 

dL 2 (u 2 ) 



d«2 



£i(ui)3Cia - L 2 (u 1 )L 1 (u 1 )X l2 'H 12 . (44) 



For simplicity, we do not write the auxiliary space dependence of the //-operator here and 
hereafter. By use of the relations 

L 1 (u)X 12 = X 12 L 2 (u), L 2 (u)X 12 = X 12 L 1 (u), (45) 

eq. (f44|) can be rewritten as 

[X u ,L 1 (u)L 2 (u)}=L 1 (u)^(u)-^(u)L 2 (u), (46) 

This is the fermionic version of the Sutherland equation, which enables us to show the 
commutativity of the Hamiltonian Ji and the transfer matrix r(u), 

FK,r(«)] = 0. (47) 

The fermionic Sutherland equation fl46|) is equivalent to the Lax equation, 



= M j+ i(u)Hj(u) - Hj(u)Mj(u). (48) 

Here Mj is given by the formula, 

M» = -ULJ\ U ){**2M + :)(,, ,.£>);}. (49) 



We remark that the formula (|4"5|) for the usual spin models is well-known. |5|] However, to 
our knowledge, the formula fl49|) for the fermion models is not known before and is proved 
here for the first time. We can conclude that the fermionic Lax equation follows from the 
Yang-Baxter relation (|26|). 

Using (|T0|) and (|25|) we can obtain the Lax operator Mj(it) for the fermion XX Z model 

as 

n ; \M 21 {u) M 22 {u)J v ; 
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where 
Mn(u) 

M 12 (u) 



1 



sin u 



1 — 



sin 2rj I sin (it + 2rj) 
i sin 4rj 



r. sinw -i | 
• I sinfit — 2??) J J_1 



sin(w + 2rf) sin(w — 2?]) 

-{in 3 -_iCj- - (1 - n^Cj-x} 



C Cj-l + — , . , 
J sm 2r? I sin(« — 2??) . 

1 — %)(1 — rij-i) — i cot (it + 2r]), 



1 



■{i(l - Uj^Cj - TljCj^}, 



sin(w — 2^) sin(it + 2r/) 

M 2 i(m) = - . 7 - - - {inj-ict + (1 - 7ij) c i-i} + o^ iK 1 - n i-i) c l + 



sin(« + 2rf) 



M 22 (w 



sm u 



sin 2r/ I sin(w + 2rf) 
i sin 4?7 
sin(it + 2rf) sm(u — 2rj) 



t 1 

J ~ J sin2?7l ' sin(it + 2?y) 



sin (it — 277) 
sin u 

1 + 



fij Tij _ 1 — icot(w + 277). 



(51) 



This result coincides with the known one |12| up to the terms proportional to the identity 
matrix. 

As another application of the fermionic Sutherland equation, we consider the boost 
operator for the fermion chain model as 



N 



3=1 



(52) 



By making use of the fermionic Sutherland equation, we can show the following relations, 

(53) 



[V,T(u)] = ^-Nstr(L N ...L 2 dLl 



du 



du 



We define I^ 1 as n-th local conserved operator and 1^ = IK Then, from (|53|), a recursion 
relation for the conserved operator is 



/(n+ i) = ^ /(n) j + ^ ome Qther terms y 



(54) 



By use of the boost operator fl5"2"|) and the recursion relation (0), we immediately obtain 
some higher conserved operators for the fermion XX Z model. Notice that the second term 



in the r.h.s. of eq. (|54]) originates from the second term in the r.h.s. of eq. (|53|), and we 
can estimate its contribution. Then we get 



1 N 

/(2) = sti T 2r^ ^ ~~ CjtCj+2) ~ 2 cos27 ?{( c i+i Cj ~ C ^' +l) 

' .7=1 



n i+2 (ct +lCj - ct Cj+ i) - (ct +2 c i+1 - ct +1 c i+a )nj 



(55) 
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2 N r 

J(3) = sh?2^ [( C i-K^ + 4 c i+3 + + c i c i+i> 



+ cos277|-3(ct +2Cj + c)c j+2 ) + 2n j+3 (c] +2 Cj + c)c j+2 ) 

+ 2n i+ i(ct +2 c,- + ctc i+2 ) + 2n J -_i(ct +2 c J - + c]c j+2 ) 

+ 2(cJ +3 c j+2 - ct +2 Cj +3 )(ct +1 c,- - ctc j+ i) - 2n 3 - + 4n i+1 n 3 - - 2n j+2 n^ 

+ 2 cos 2 2»7|(ct +1 c^ + ctcj+i) + 2n j+3 (c] +2 c j+1 + cj^c,-^)^ 

- n i+2 (ct +1 c,- + c]c i+ i) - nj(c] +2 Cj+i + c] +l c j+2 ) } . (56) 



We remark that these results are consistent with the known ones. 11, 12, 19 



5 Concluding Remarks 

In this paper, we have studied the quantum inverse scattering method (QISM) for the 
XX Z fermion model. We have introduced a fermionic /^-operator for the fermion chain 
model. It consists of the fermion operators and satisfies the Yang-Baxter equation. The 
i?-operator intertwines the fermionic L-operator in a different way from the usual one. 
The auxiliary space and the quantum space of the fermionic L-operator are exchanged. 
We have shown that the new relation plays a complementary role to the usual graded 
Yang-Baxter relation. In particular, we can derive the fermionic Sutherland equation from 
the relation. It is shown that the fermionic Lax equation is equivalent to the fermionic 
Sutherland equation. The Sutherland equation also gives a mathematical foundation of 
the boost operator approach for the fermion chain model. 

We have also defined a fermionic monodromy operator by a product of the fermionic 
i?-operators. The transfer operator, which is an analog of the transfer matrix, can be con- 
structed by taking the expectation values in the auxiliary space, Str, of the monodromy 
operator. The transfer operator constitutes a commuting family, which proves the integra- 
bility of the fermion chain model. The diagonalization of the fermionic transfer operator 
will be discussed in a separate paper. 

Our approach is applicable to other integrable fermion models. Among them, we shall 
consider the Hubbard model in subsequent papers. 
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